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We show that patterned defects can be used to disrupt the sub-lattice symmetry of graphene so 
as to open up a band gap. This way of modifying graphene's electronic structure does not rely on 
external agencies, the addition of new elements or special boundaries. The method is used to predict 
a planar, low energy, graphene allotrope with a band gap of 1.2 eV. This defect engineering also 
allows semiconducting ribbons of carbon to be fabricated within graphene. Linear arrangements of 
defects lead to naturally embedded ribbons of the semiconducting material in graphene, offering the 
prospect of two-dimensional circuit logic composed entirely of carbon. 



Graphene-based electronics promise to overcome 
limitations associated with silicon technologies [THJ] 
and open previously unavailable new applications [5]. 
The improved performance derives from unparalleled 
strength [6] combined with high electron mobility [3j [7] . 
Graphene is a semi-metal, though, and therefore requires 
modification to open up the band gap desirable for logic 
circuits. This is simple to do in principle, as demon- 
strated in Figure [T] where small rotations of a bond, peri- 
odically repeated through the lattice, collapses the Dirac 
cones at the graphene K point. Continued rotation of 
the bond will close the band gap elsewhere resulting in 
a conductor. To exploit this sensitivity, an agency must 
be found for tuning the distortion of lattice symmetry. 
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FIG. 1: The twisting of a bond within a periodic cell (a) 
disrupts lattice symmetry, causing band gap opening (b,c) 
and closing (d). The Fermi energies have been aligned to 
more clearly see the effect of bond rotation. 



A number of approaches have been considered for dis- 



torting graphene so that it exhibits a band gap. All of 
these rely on either the action of an external agency or 
a confining boundary to disrupt the delicate sub-lattice 
symmetry. These include the use of strongly interacting 
substrates [8HTD] or the application of electric fields in 
graphene bi-layers [8] [TTJ [12] . Quantum confinement ef- 
fects can also be exploited as in carbon nano-ribbons [T3l - 
[16], or carbon nano- meshes [17]. The graphene can also 
be chemically modified. For instance, hydrogen can be 
introduced so as to change the bonding from sp-2 to sp-3 
to create the corrugated insulator graphane [18]. Simi- 
larly, a graphene oxide insulator can be created via bond- 
ing with hydroxyl groups [19] . The external agency could 
also be in the form of neighboring regions of foreign ma- 
terial as in composite two-dimensional structures that 
harbor domains of graphene [20] . 

Here we explore an alternative approach wherein the 
symmetry of the graphene is modified using patterned de- 
fects. Previous defect investigations have predicted only 
metallic allotropes of graphene [21-24 . Interestingly, 
all of these conducting allotropes can be constructed 
by patterning two types of defects [23], Stone-Thrower- 
Wales (STW) [25] [21] and Inverse Stone-Thrower- Wales 
(ISTW) [271 US defects. The consistent metallic na- 
ture of these graphene allotropes, though, suggests that 
their local lattice distortions are not appropriate for del- 
icately teasing apart the K point Dirac cones without 
pinching off the band gap elsewhere. This motivated the 
search for a new defect which might produce semicon- 
ducting graphene, and we have discovered that the di- 
vacancy (DV) defect is well-suited for the task. A pat- 
terned combination of DV and STW defects results in 
a planar structure with a band gap, a semi-conducting 
form of graphene. This stable, single-layer, all carbon 
semiconductor, named Octite SC, is shown in Figure [2] 
We recently became aware of a complimentary approach 
in which pentagon and heptagon defect were coupled 
with surface corrugation to introduce a band gap in 
graphene [29] . 

As a starting point, Density Functional Theory (DFT) 
was employed to estimate the ground state structure 
and predict electronic character. A real-space numer- 
ical atomic orbital code [30 was used to initially re- 
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FIG. 2: A computationally created semiconducting allotrope 
of graphene resulting from patterned defects. The planar 
structure is only 313 meV/atom above graphene, comparable 
to 233 meV/atom for Haeckelite H567, a metallic graphene 
allotrope previously theorized [22] . 



lax the structures until the energy change was less than 
2.7 x 10 -4 eV. The Perdew-Wang generalized gradient 
approximation accounted for electron exchange and cor- 
relation energy [31] . This computational setting has been 
shown to accurately predict defect geometries and ener- 
gies associated with graphene structures [23l [271 12H] • 

The new allotrope has a primitive cell with 28 atoms 
composed in pentagons, hexagons, and octagons. It has 
a planar density of 0.364 atoms per A 2 , comparable to 
0.380 atoms per A 2 for graphene. The central geomet- 
rical feature is an octagon completely surrounded by 
hexagons. The cell geometry is square with p4/MMM 
symmetry. A Hessian linear vibrational analysis showed 
the structure to be a local energy minimum. Room tem- 
perature quantum molecular dynamics simulations were 
also used to confirm its stability. 

As a final check on the ground state structure, the 
cell geometry and linear stability were verified using a 
second DFT implementation that employs a plane-wave 
basis set [32j [33]. A Projector Augmented Wave ap- 
proach [34] was used within a Generalized Gradient Ap- 
proximation [35]. The wave function energy cutoff was 
set at 400 eV, and a conjugate gradient method employ- 
ing a 4x4x2 Monkhorst-Pack grid resulted in the same 
ground state structure but with a lattice constant of 8.674 

A. 

The plane-wave DFT code was subsequently used to 
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FIG. 3: Band structure (BS) and density of states (DOS) for 
a semiconducting graphene allotrope. Both DFT and GoWo 
predictions are shown. The DOS plot at right is from the 
DFT calculation. 



predict both the band structure (BS) and Density of 
States (DOS). A 24x24x1, T-centered k-space grid was 
used to calculate the electronic properties. The BS and 
DOS diagrams are shown in Figure [3J The material pos- 
sesses a 0.2 eV direct band gap at the T point. As 
in graphene, all carbon atoms bond with three nearest 
neighbors, but the symmetry of the in-plane sp-2 char- 
acter has been disrupted by the varying bond angles and 
bond lengths. 

DFT typically underestimates the band gap of periodic 
structures [36], so the 0.2 eV gap should be viewed as 
a lower bound. Many-body perturbation theory, within 
the GqWq approximation [37] . was therefore used to pro- 
vide a more accurate estimate of the electronic struc- 
ture. As expected, the inclusion of quasi-particle effects 
widens the band gap, providing an estimate for the di- 
rect, T point gap of 1.1 eV. The GqWq band structure 
is shown with its DFT counterpart in Figure [3] Three 
iterations towards self-consistency (GWo approximation) 
opened the band gap slightly to 1.2 eV. 

Having elucidated the lattice and electronic structure 
of this new material, we now examine its defect archi- 
tecture. Such decompositions are not intended as pre- 
scriptive synthesis steps and are, rather, a way of un- 
derstanding the properties of graphene allotropes from 
the broader perspective of defects from which they are 
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FIG. 4: Patterned defect composition of Octite SC from 
graphene. (a) Graphene with a di- vacancy, (b) A STW defect 
is introduced by rotating the gray (green) bond, (c) A second 
STW defect is introduced, (d) Patterning of the patch re- 
sults in a second octagonal defect surrounded by alternating 
hexagons and pentagons, shown in gray (red). 



composed. All Haeckelite allotropes can be constructed 
from graphene with templates of either STW defects or 
a combination of STW and ISTW defects [23]. Like- 
wise dimerites, as the name suggests, amount to graphene 
on which a regular grid of ISTW defects have been ap- 
plied [23 . The new semiconducting allotrope, however, 
requires the addition of di-vacancies (DV) to the defect 
alphabet. By patterned placement of DVs and STW de- 
fects, Octite SC can be engineered from graphene. 

We schematically illustrate the defect composition in 
Figure 4 Initially a DV defect is introduced as shown in 
Figure 4 'a). Two bond rotation defects are subsequently 
imposed (frames b and c). The result is then replicated 
so as to create an array of octagonal defects. Octagons 
surrounded by hexagons result from DV and STW de- 
fects while smaller octagons surrounded by alternating 
pentagons and hexagons derive from STW defects (frame 

d). 

The formation energy associated with each defect 
shown in Figure [4] was calculated within a periodic 200 
atom supercell as well as on a passivated, 206-atom flex- 
ible graphene flake. The formation energy of the DV 
was found to be 7.6 eV in the periodic supercell, in good 
agreement with previous periodic DFT estimates of 7.5 
eV [38] and 8.7 eV [39]. The flexible boundaries of the 
graphene flake allowed the formation energy to be sub- 
stantially lower at 6.6 eV. We calculate that the subse- 
quent bond rotation defects on the right and left of the 



DV result in sequential, net increases of 4.7 eV and 5.1 
eV for the periodic supercell or 5.0 eV and 5.1 eV for the 
graphene flake. 

DV defects, central to the structure of Octite SC, can 
be precisely created via irradiation [40]. Recently, an 
electron irradiation beam focused to 1 A was shown to 
be capable of precisely creating vacancies in CNTs [41]. 
The stability of DVs has also been considered theoret- 
ically within a CNT setting, where molecular dynam- 
ics annealing predicted that they will remain stable in 
the 5-8-5 configuration up to at least 2700 K [42^ and 
3000 K [43]. Tight-binding molecular dynamics studies 
at higher temperatures concluded that DVs were stable 
out to 3100K for at least the 90 ps computational exper- 
iment [44]. A subsequent temperature increase to 3800 
K was required to cause this 5-8-5 defect to relax to the 
slightly lower energy 555-777 defect with an estimated 
barrier of 5.2 eV [44 . Our own transition state analysis 
estimated a similar barrier of 3.7 eV for relaxation from 
a 6-8-6 defect to a 555-777 defect. These results sug- 
gest that the DV defect, central to Octite SC, should be 
stable for the operating conditions expected of standard 
electronic circuitry. 

STW defects have been extensively studied in both 
graphene [25] [45] and CNTs [38] [46]. STW defects could 
be created or reversed by scanning tunneling microscopy 
(STM) [47] or atomic force microscopy (AFM) [45] . 
Both DV and STW defects have been experimentally ob- 
served @oi sh [49]. 

Semiconducting ribbons can be built into graphene as 
shown in Figure [5] Our DFT calculations predict that 
these ribbons remain planar within the graphene, and 
preliminary results suggest that the semiconducting char- 
acter is locally preserved. Analogous conducting struc- 
tures can also be designed, and our prediction of a metal- 
lic defect ribbon has already been experimentally veri- 
fied [24] [50]. Semiconducting and conducting ribbons, 
embedded within graphene, could be used together as 
the basis for electronic logic circuits. 

In summary, we have found that patterned defects in- 
volving DVs yield a planar, single- atomic layer material 
composed entirely of carbon that is an intrinsic semi- 
conductor; no external agencies or special boundary con- 
ditions need to be imposed in order to open up a band 
gap. In comparison with epitaxial paradigms, this is par- 
ticularly attractive since suspended graphene possesses 
a mobility ten times that of graphene fabricated on a 
substrate [3] . Future electronics designs can exploit local 
and extended defects to build free standing, single-atomic 
layer logic circuits based on the local symmetry proper- 
ties of carbon bonds. 
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FIG. 5: A ribbon of Octite SC built within graphene using 
DV and STW defects. 
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